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Problem  and  Results:  Technical  Report 

Contract  DAAG29-83-K-0025 

Opt iaal  Design  of  Fibered  S tructures 

Professor  Gilbert  Strang 
Department  of  Mathematics 
Massachusetts  Institute  of  Technology 


The  classical  variational  problems  involve  the 
minimization  of  convex  energies  subject  to  boundary  conditions. 
They  lead  to  elliptic  differential  equations  and  the  solutions 
are  smooth  functions  away  from  singularities  in  the  data  or  the 
domain.  A  typical  equation  is  div  (c(x)  grad  u)  =  f,  with  a 
coefficient  c  that  reflects  the  material  properties. 

In  contrast,  the  problem  of  optimal  design  is  to  choose 
that  coefficient  c.  The  cost  function  may  still  be  quadratic 

in  u,  or  in  plasticity  problems  it  may  involve  the  L**  norm 

of  c  grad  u,  but  the  choice  between  different  materials  or 

different  shapes  is  determined  by  the  unknown  function  c:  it 

is  a  problem  in  which  the  coefficient  is  the  control  variable. 

In  such  problems  it  is  common  to  find  that  a  conventional 

solution  does  not  exist.  The  minimum  cost  is  approached  by  a 

sequence  c  that  is  more  and  more  oscillatory.  That 
n 

corresponds  to  a  series  of  designs  in  which  materials  are 
combined  with  an  increasingly  complicated  microstructure. 

Mathematically  the  underlying  difficulty  is  that  the 
variational  problem  is  not  convex.  In  optimal  design  the  cost 
may  be  zero  in  the  absence  of  material  and  jump  to  unity 
whenever  the  stress  is  nonzero  and  material  is  required.  The 
minimizing  sequence  of  coefficients  makes  that  jump  more  and 
more  frequently,  and  there  is  no  convergence  in  norm  to  a  limit 
c  (x).  However  there  is  weak  convergence.  The  quantity  of 
material  in  each  region  does  approach  a  limit,  and  the 
orientation  can  be  defined.  The  limit  can  be  described  as  a 


materials . 


The  mathematical  theory  is  expressed  by  a  relaxation  of 
the  original  nonconvex  problem.  The  problem  is  replaced  by  a 
well-posed  (lower  semicont inuous )  problem  which  has  the  same 
minimum  cost,  but  that  minimum  value  is  actually  achieved  (by  a 
composite  material).  The  solution  of  the  relaxed  problem  is 
the  weak  limit  of  minimizing  sequences  in  the  original  problem. 
Our  joint  work  with  Robert  Kohn,  referenced  in  the  bibliography 
below,  carried  out  an  explicit  relaxation  and  solution  of 
problems  in  optimal  design.  It  is  comparatively 
straightforward  in  the  case  of  scalar  problems,  but  it  becomes 
more  complicated  when  the  unknown  u  is  a  vector  function. 
Elasticity  and  continuum  mechanics  produce  the  vector  case,  and 
much  of  our  research  was  concentrated  on  the  relaxation  of 
functionals  JF(grad  u)dx,  in  which  grad  u  is  the  Jacobian 
matrix  and  ouasiconvexitv  is  the  key  property  to  achieve  for 
F. 

The  long  paper  [1]  summarizes  our  work  on  the  interaction 
of  homogenization,  relaxation,  and  the  mathematical  theory  of 
composite  materials.  In  this  report  I  will  describe  one 
specific  problem  solved  separately  in  [4].  It  was  a  central 
problem  in  Tartar’s  development  of  the  method  of  compensated 
compactness,  and  it  is  an  optimization  problem  for  the 
effective  properties  of  composites:  How  should  conducting 
material  be  placed  in  a  unit  square  so  as  to  carry  prescribed 
currents  in  the  north-south  and  east-west  directions,  with  a 
minimum  area  of  material?  We  formulated  this  question  as  a 
nonconvex  variational  problem: 

2 

Minimize  the  area  of  the  subset  S  c  [0,1]  subject  to 
2  2 

//|vu|  dxdy  =  0,  //Jvvj  dxdy  =  R,  vu  =  vv  =  0  outside  S, 
with  boundary  conditions  u  =  0  on  the  left  side,  u  =  1  on 
the  right,  v  =  0  on  the  bottom,  and  v  =  1  on  the  top.  The 


streaa  functions  are  u  and  v  (a  vector  unknown),  the 
effective  resistances  of  the  square  are  Q  and  R,  and  the 
problea  is  to  achieve  those  resistances  with  a  ainiaua  of 
aaterial.  The  area  ( Q+R-2 ) / (QR-1 )  can  be  achieved  by  a 
suitable  aicrost ructure  and  the  problem  is  to  prove  that  this 
area  is  minimal. 

The  nonconvexity  enters  in  the  expression  /Jig  ^x<*y  f°r 
the  area;  the  characteristic  function  equals  one  in  S  and 
zero  in  the  complement.  The  difficulty  is  to  compute  the 
" quas iconvex i f icat ion"  of  the  functional  F(vu,vv)  = 

2  2 

lg  +  a|vu|  +  |  w  j  ,  in  which  the  Lagrange  multipliers  A  and 

/j  enter  from  the  constraints.  Xohn  and  I  determined  the 
largest  function  PF  below  F  that  is  polyconvex  (convex  in 
vu.vv,  and  det  J).  We  also  established  that  this  function  is 
convex  on  lines  connecting  Jacobian  matrices  J  =  [vu  vv]  that 
differ  by  rank  one.  Therefore  PF  is  the  correct 
quas iconvex i f i cat  ion  of  F,  and  the  relaxed  problem  requires 
the  minimization  of  //PF< vu , vv ) dxdy .  That  calculation  led 
back  to  ( Q+R-2 ) / ( QR- 1 )  as  a  lower  bound  on  the  area,  and 
therefore  the  ainiaua. 
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